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Abstract 

^: 

1^ . We study the effects of including one-loop radiative corrections in a non-super- 

symmetric hybrid inflationary model. These corrections can arise from Yukawa cou- 
plings between the inflaton and right-handed neutrinos, and induce a maximum in 
^ . the potential which admits hilltop-type solutions in addition to the standard hybrid 

lO \ solutions. We obtain a red-tilted spectral index Ug, consistent with WMAP5 data, for 

' sub-Planckian values of the field. This is in contrast to the tree level hybrid analysis, 

Ti^lj- ■ in which a red-tilted spectrum is achieved only for trans-Planckian values of the field. 

Successful reheating is obtained at the end of the inflationary phase via conversion of 
' the inflaton and waterfall fields into right-handed neutrinos, whose subsequent decay 

■ can explain the observed baryon asymmetry via leptogenesis. 

It has been shown in the recent WMAP five year analysis (WMAP5) that tree level 
^ ■ chaotic inflation driven by a quartic potential is excluded at the 99% confidence level [1]. 

■ The tensor to scalar ratio r in such a model turns out to be excessively large. However, this 
conclusion can be significantly altered provided one takes into account radiative corrections, 
especially those generated by any Yukawa couplings that may be present between the 
inflaton and fermion fields in the model. Such Yukawa couplings can be expected on 
general grounds, particularly in those inflation models containing right-handed neutrinos 
in which the observed baryon asymmetry is explained via type I leptogenesis [2]. It has 
been shown that radiatively improved quartic potential models with plausible values of the 
Yukawa couplings can yield values for r and the scalar spectral index Ug that lie well inside 
the 2-a bounds provided by WMAP5 [3j. 

In this paper, following Ref . [3] , we carry out a similar analysis for a hybrid inflationary 
(HI) potential [4]. In addition to providing a mechanism for generating the primordial 
baryon asymmetry, the extra couplings of the scalar fields in the model (especially the 
inflaton) to right-handed neutrinos also play an important role in how the scalars decay, 
and thus in the reheating of the universe. We will see that, as in the quartic case, the 



1 



radiatively corrected model explores an expanded region of parameter space, and can result 
in better agreement with experimental data. 

At its tree level, the HI potential can be written as [5] 



where M, m are mass parameters and k, A are dimensionless. The global minima of the 
potential lie at ((x); (0)) = (±2M, 0). The effective mass squared of the field x iii the 
X = direction is = — k^M^ + A^0^/2. Thus, for (p > (pc = the only minimum of 

the potential V{x, (p) lies at x = 0. In this region the HI potential takes the form 

v{<i>) = v, + = K) [l + , (2) 

where = J!^y^ , and Vq = k^M^ is the constant vacuum energy term. The second term 
in Eq. ([2]) provides a non-zero slope in the otherwise flat potential, and the system can 
inflate as it rolls down the x = valley. Upon reaching cj) = (p^, the minimum in the x 
direction becomes a maximum, and inflation ends abruptly as the system rapidly falls into 
the global minimum. This scenario is termed "hybrid" because the vacuum energy density 
Vq is provided by the waterfall field x, while is the slowly rolling infiaton field. 

Before considering radiative corrections, it is appropriate to discuss the predictions of 
the tree level hybrid inflationary (TLHI) model in comparison with WMAP5. The slow- roll 
parameters for the TLHI potential are given as 
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where r]o = ri{(p = 0) = — "^-P — 2.4 x 10^^ GeV is the reduced Planck mass. The 
number of e-foldings after the comoving scale Iq = 2TT/ko has crossed the horizon is given 
by 

^^0 = ^ rYT^)#=- r 175^ M0 = ^(0o-0c + ln 0^/0^), (5) 

where 0o is the value of the field when the scale corresponding to exits the horizon, and 
0c = sJ^'^Q^'^ (or 0c = {(pc/mp)), with K ~ A ~ 10~^, is the value of the field at 
the end of inflation. In Eq. ([5]), we may eliminate r/o iii favor of Vq using the curvature 
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Figure 1: r ys Ug for the tree level hybrid inflationary potential, shown together with the 
WMAP5+BA0+SN contours (68% and 95% confidence levels) [1]. Here the small (large) 
dots correspond to Nq = 50 (A^o = 60). 
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In our calculations, we will use the value AjilkQ) = 4.91 x 10^^ obtained by the recent 
WMAP5 analysis for ko = 0.002 Mpc~^ [1]. To leading order, the spectral index and 
the tensor-to-scalar ratio r are given by 

- 1/2 

Ug — 1 — 6e + 2?7 = 1 — 4 ?7o 
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Using the above equations, we show in Fig. [T] the predictions of the TLHI model along 
with the WMAP5 1-a and 2-a bounds, and specify the range in which 50 < A^o < 60. 
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The flat potential regime (0o ^ 1) lies outside the 2-cr bound, corresponding to a blue 
spectral index. A red-tilted spectral index is obtained for 0o > l/"\/2 — 0.71, and for even 
larger values of 0o this result reduces to the quadratic potential prediction r ~ 4(1 — Ug). 
It is interesting to study the value of the field (p compared to the reduced Planck mass 
mp (i.e., (j)p = ^) for which the sprectral index is red-tilted. With 0o > l/"\/2 and the 
slow-roll approximation (r], e) <^ 1 (which implies ?7o <^ 1), the blue spectrum can only be 

avoided for trans-Planckian values of the field (hp = \ — (t>n > \ — ^ 1- For a more exact 

treatment obtaining a red spectrum in the TLHI model, see the recent paper in Ref. [B]. 
With these tree level results in mind, we would now like to examine the situation after 
including radiative corrections in the HI potential. 
Consider the following Lagrangian density 
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- \yB(t^BNN - WbXbNN - ^m^viViV, (10) 



where the subscript 'B' denotes bare quantities. To keep the discussion as simple as 
possible, we have introduced a single Yukawa coupling involving and each of and 
X- In a more realistic scenario, successful leptogenesis requires at least two right-handed 
neutrinos. Since x = during inflation, the interaction between and x has no effect in 
this regime. After inflation, however, A^ acquires a contribution M^v — Y {x) to its mass in 
addition to the bare mass uin. Also, oscillations of x affect the way in which the system 
reheats after inflation, as we will see later when we discuss the reheating phase. 

The inflationary potential including one-loop corrections, in terms of renormalized 
quantities, is given by 
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where Vioop is the one-loop correction to the tree level potential. In the x 
Vioop can be written as [7J 
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During inflation, (j) is always larger than <^c, therefore for ycjyc ^ {mj\f,m) and y > the 
above potential reduces to 
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In order to ensure that the log factor is always positive during inflation, it is convenient to 
set the renormalization scale at /i = yipc- Then, the complete radiatively-corrected hybrid 
inflationary (RCHI) potential in the x = direction reduces to the form 

.2 J,2 
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proximate A^ to be independent of 
by 
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In the following calculations, we will ap- 
The slow-roll parameters in this case are given 
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Using Eq. ([6]) and solving for rjo as before, we find 
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The number of e-foldings in the RCHI model can be calculated as 
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In order to ensure that the potential remains bounded during inflation, we take 0c < 0o < 
^ where (j)M is field value at the maximum of the potential. This maximum 
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introduces hilltop-type solutions [8J, to which we now turn our discussion. 

In its approximate form, the potential in Eq. (IT^ has previously been analyzed only 
for sub-Planckian hilltop-type solutions [9] , for which inflation begins in a region where the 
potential is concave downward. In general, the RCHI model can lead to both hilltop and 
non-hilltop solutions (see Ref. [3] for a similar analysis in the case of radiatively-corrected 
quadratic and quartic potentials). In order to study both types of solutions, it will be 

convenient to define the quantities / = 4>q/4>Mi fc = 4'c/4>m and fi = \ 2 A^. As we will 



see, for the same values of A^ and Vq, we can always distinguish two separate branches 
of solutions, one with large / values and another with small / values. To facilitate our 
discussion, we can rewrite the number of e-foldings in the form 
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Figure 2: Nq and logio [M vs / for /i = 0.1, 1 and 10 with V^^^ = (lO^^-^^ - lO^^-^^) GeV, 
(^^Qi6.08 _ 1016-22) Qgv and (lO^^-^i - 10^^-^^) GeV respectively. We obtain two solutions 
(ns,r) for each value of and Vq. 



In order to discuss the two branches qualitatively, we note that the factor inside paren- 
theses in A^o is slowly varying over the entire range of /. Therefore, it is sufficient to 
consider only the contribution of the t]q^ factor in Nq. For large values of /i, r]Q^ reduces 
to the form /^(l — /^)^, with a maximum at / = 1/a/3 ~ 0.58. As we move away from 
this maximum, we obtain the same value of Nq for each of two values of /, one smaller and 
one larger than These large-/ and small-/ branches can be seen in Fig. [2l As /i 

is lowered, the maximum of A^o moves toward smaller values of /, asymptotically tending 
toward / = fi/V2. 

Let us discuss in greater detail the large /i limit, which corresponds to the "flat potential 
regime" (i.e. 0o = ///i -C 1). In this limit, the number of e-foldings A'o reduces to 



i^J^.]f- (21) 



with 



l2yTi^/\l \ mp ) \Pil-f 

An upper bound on Vq can be found by considering the point at which the two branches 
meet for a given value of /i. Owing to the inverse dependence of A^o on K)? the number 
of e-foldings at the maximum is shifted upward as Vq is decreased away from its limiting 
value, and the two branches diverge from one another in the range corresponding to realistic 
values of A'q. The basic condition f > fc (or A'o > 0) leads to the bound Vq^'^ < ^^^j^ 
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Figure 3: r vs n^, for radiatively-corrected hybrid inflation, with / = 0.4 and / = 0.65. 
Here, small (large) dots correspond to A^o = 50 {Nq = 60). Compared to Fig. [H smaller 
values of <^o fell within the WMAP5 bounds. In addition, the allowed and disallowed 
regions are exchanged in going from a non-hilltop (/ < /j) to a hilltop (/ > /j) solution. 



GeV. In order to obtain a realistic number of e-foldings, Vq must obey a somewhat more 
stringent bound. 

As mentioned earlier, hilltop solutions are deflned as having a concave downward cur- 
vature of the potential when inflation begins. Therefore, only those solutions which satisfy 
/ > /j = l/-\/3 ~ 0.58 will be called hilltop solutions, where is the value of / at the 
point of inflection. 

The spectral index Ug and the tensor-to-scalar ratio r for the RCHI model are given as 



l-2(l + 3A^)02 + 5A^04_g^2^6 

= l + 2r]o^ z ^-TTl (22) 
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r = ^ ^^=m-ns)^ — ^-^ . (23) 

The spectral index becomes unity when 0o acquires the value 




(24) 



which reduces to the tree level result (pQ = in the f —>■ limit. For 0o values larger 
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Figure 4: r vs Ug for radiatively-corrected hybrid inflation, with (pM = 100, 1 and 0.1 for 
/ = 0.4, 0.5, /j, 0.7 and 0.8. Large values of 0m result in appreciable values of r for a 
realistic number of e-foldings. 



(smaller) than in Eq. (p^ . we obtain a red-tilted (blue-tilted) spectral index if / < /«. For 
f > fi, r] becomes negative and we always have a red-tilted spectral index. 

For small values of /, Eq. ( l23l) reduces to the tree level result r ~ 4(1 — Ug) -J^^ , ■ 

'Po~ / 

Similarly, for large values of (po we find r ~ 4(1 — Us) , 5 2 3 4V ^'^i'^^ reduces to the 

quadratic prediction r ~ 4(1 — n^) in the small / limit (see Figs. OandHj). Thus we see 
that in order for the radiative corrections to produce a reasonably large contribution, 0o 
should be chosen close to ipM- This choice in turn shifts the "fiat potential regime" toward 
the inside of the 2-0" bound of the WMAP+BAO-I-SN data, as can be seen by comparing 
Figs. [1] and [31 A larger choice of / leads to a greater portion of the "fiat" region inside 
the bounds. Additionally, Fig. [3] shows that for f < fi the allowed regime lies above the 
00 = 100 line, whereas for f > fi the allowed regime lies below this line. As noted earlier, 
for f > fi, f] becomes negative and both e and 77 tend to drive Ug {— 1 — Q e + 2 rj) below 1 
as ?7o increases. In contrast, for f < fi, rj is always positive and competes with e, allowing 
Us values above or somewhat below 1. 

In our discussion so far, we have suppressed the dependence of in order to obtain 
analytically tractable expressions. Next, we employ numerical calculations to examine the 
RCHI model in such a way that this dependence can be taken into account. In these 
calculations, we use the next-to-leading order expressions for n^, r, and A-ji for added 
precision [12]. Having already explored the relative behavior of the quadratic and vacuum 
terms, we will use the form of the potential Eq. fll4p written in terms of rather than 
in order to more directly probe the parameters and Vq. 

As discussed earlier, inflation ends via a waterfall effect at the critical value 0c. It is 
worth noting that inflation can end due to the breakdown of the slow-roll approximation, 
before reaches this value; however, in our calculations below, |?7| < 10~^ and e < 10~^ 
when evaluated at 0c, and so the slow-roll relations are still valid when the waterfall is 
induced. 
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In terms of the reheat temperature T^, the number of observable e-foldings may be 
written as [10] 



where Tr 
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y/Trrip [TT], and F is some appropriate decay width. 

Two sources of primordial production contribute to reheating, resulting from the oscil- 
lations of and X about their respective minima. The decay rates of these two processes 
are given by 



r 



X^NN 



(j>-*N N 



y 



Stt ' ^ '''' Stt 
where Y ~ Mn / (x) ; and the scalar field masses are given as 



(26) 



V2kM, = sjm? + 2{\Mf. (27) 
To see how these two contributions compare, consider the ratio of their decay rates 
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<t>~>N N 
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(28) 



To simplify the analysis, we will take A ~ /t ~ ?/ ~ F. Then, for m <C M, we obtain 
~ and the ratio of the decay rates in Eq. fl28|) is of order unity. Indeed, under 
these assumptions, it turns out that all relevant decay widths are of the same order, and 
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Vr (GeV) 


.4(10-1'*) 


M (GeV) 


m (GeV) 


fc 


n,. 


r(10-'') 


777^(10"GeV) 


T, (1012 GeV) 


Na 


<1>P = 0.25 


0.65 


3.696 X 10^" 


1.164 


5.744 X 10"' 


9.04 X 10" 


0.1735 


0.9864 


0.149 


1.68 


1.45 


56.30 


0.70 


3.584 X lO'-'^ 


1.282 


5.503 X 10"' 


8.74 X 10" 


0.1789 


0.9717 


0.132 


1.65 


1.47 


56.29 


0.75 


3.434 X 10'-^ 


1.336 


5.245 X 10"' 


8.29 X 10" 


0.1828 


0.9536 


0.111 


1.59 


1.46 


56.25 


0.80 


3.235 X 10'-' 


1.296 


4.960 X 10"' 


7.63 X 10" 


0.1843 


0.9312 


0.0870 


1.49 


1.41 


56.19 


4>P = 1 


0.60 


8.873 X 10'-' 


2.134 


1.278 X 10'' 


6.13 X 10'^ 


0.08907 


0.9930 


5.02 


4.36 


2.72 


57.10 


0.65 


8.603 X 10'-^ 


2.333 


1.225 X 10'^ 


5.87 X 10'^ 


0.09252 


0.9768 


4.43 


4.27 


2.75 


57.08 


0.70 


8.262 X 10^-^ 


2.430 


1.171 X 10'^ 


5.53 X 10'2 


0.09521 


0.9575 


3.77 


4.12 


2.73 


57.04 


0.75 


7.836 X 


2.387 


1.113 X 10'^ 


5.10 X 10'2 


0.09696 


0.9346 


3.04 


3.90 


2.65 


56.99 


dp = 2.5 


0.55 


1.482 X 10'" 


2.280 


2.116 X 10" 


1.88 X 10" 


0.05409 


0.9921 


43.1 


7.34 


3.59 


57.55 


0.60 


1.443 X lO^** 


2.549 


2.032 X 10'^ 


1.81 X 10" 


0.05665 


0.9783 


38.7 


7.24 


3.67 


57.53 


0.65 


1.394 X lO^" 


2.728 


1.947 X 10'^ 


1.72 X 10" 


0.05880 


0.9619 


33.7 


7.06 


3.68 


57.50 


0.70 


1.334 X 10'" 


2.779 


1.859 X 10'^ 


1.60 X 10" 


0.06046 


0.9427 


28.2 


6.77 


3.62 


57.48 



Table 1: Predicted values of various inflationary parameters in the radiatively-corrected 
hybrid model. Here we show only those values falling within the WMAP5+BA0-I-SN 2- 
cr bounds (see Fig. E]) for each of the three cases 0p = ^ = 0.25, 1, and 2.5, which 
correspond primarily to hilltop solutions (/ > 0.6). These values correspond to the choice 

0,~2\/oVv^in Eq.dSO]). 
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Figure 5: logigA vs. / and Us for the cases 0p = = 0.25, 1, and 2.5. Two solutfons 
of each of /, Ug exist for a single value of A (see Fig. [2]). The region falling within the 
WMAP5 bounds mainly leads to two hilltop solutions. 



so we may approximate the total decay width appearing in as F ~ F^^tv n- Then, using 
Wreh = for matter dominant reheating and taking ~ 106, the reheat temperature 
becomes ~ 0.035 j/^/m^ mp. 

To proceed further, we may eliminate in favor of the field value at the local max- 
imum (pM induced by the radiative correction term. Assuming this maximum is the only 
extremum other than the minimum at the origin, we can write 



Using the approximations above, we may simplify the expression for the critical value of 
the inflaton field: 



Note that this expression depends only on Vq and A (via y) . 

We are interested in comparing sub-Planckian and trans-Planckian inflation. We thus 
consider three values of the inflaton at the start of inflation: 0p = 0.25, (pp = 1 and 
(f)p = 2.5. In each of these cases, if the ratio / is fixed, (pM is known and the inflationary 
potential V{(j)) is specified in terms of Vq and A. 

To perform our numerical calculations, we fix a value of / and scan over values of Vq 
until the number of e-foldings Nq given by an integral similar to Eq. ([5]) matches its value 
as given by thermal considerations, Eq. fl25|) . For each (0o, /, Vq), the value of A can then 
be calculated by setting the curvature perturbation equal to its WMAP5 value. The results 
of these calculations are displayed in Table [H 

The RCHI model yields values of the reheat temperature on the order of ~ 10^^ GeV, 
which is substantially larger than the range allowed by supersymmetric models of infiation. 




(29) 




(30) 
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Figure 6: log^g'" ^s. Ug for the cases 0p = ^ = 0.25, 1, and 2.5, shown together with 
the WMAP5+BA0+SN contours. The smaller black dots denote / values increasing from 
right to left in increments of 0.05, and the larger colored dots represent inflation starting at 
the inflection point of the potential, with /j ^ 0.60. With larger values of 0p, the WMAP5 
limit on r can be saturated. 



Tr ~ 10^-10^ GeV [13]. This range is formulated based on the gravitino constraint, which 
is peculiar to susy models and does not pose a danger to the non-susy model that we 
currently consider. Table [1] also shows that, in order to obtain a spectral index near the 
WMAP5 central value, the RCHI model predicts Vq^'^ ~ 10^^-10^^ GeV. A similar range 
of preferred Vq^^ values has recently been realized in other models such as susy hybrid 
inflation [H] and non-susy Coleman- Weinberg and Higgs inflation [15j. In the non-susy 
models, this range was also seen to be associated with proton decay, predicting a lifetime 
of 10^^-10^^ years. In addition, such models predict a reheat temperature 4-6 orders of 
magnitude lower than the RCHI model, due in part to a Yukawa coupling on the order of 
10~^. In contrast, the Yukawa coupling involving (p in the RCHI model is linked to the 
coefficient A ~ 10~^^, implying y ~ 10"^. 

Fig. [5] shows the behavior of A with respect to both / and n<j in each of our three 
cases. Each value of A corresponds to two solutions of / and Ug, as expected based on 
our previous calculations. Given our current choice of fixed parameters, some values of A 
result in one hilltop solution and one non-hilltop solution, while others result in two hilltop 
solutions. In the numerical case, the value of fi deviates somewhat from its approximate 
analytical value l/y/S, hovering close to 0.60 for all three cases. 

The results of our numerical calculations for the RCHI model are compared to WMAP5 
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in Fig. [6l As / increases along the curves of constant 0p, Ug decreases to values less than 
unity. If / is increased further, these curves pass into the region favored by WMAP5, 
span the full 2-a range, and pass outside the bounds again. For the values oi (j)p which we 
have considered, we obtain a favored range of / corresponding to 0.55 < / < 0.80. Stated 
differently, the inclusion of suitable radiative correction effects in the HI scenario results 
in the model becoming more favored by the WMAP5 data. 

In the range 0.4 < / < 0.9 that we have explored, each choice of produces a distinct 
prediction for the size of r. The case (f)p = 2.5, which admits trans-Planckian values of the 
inflaton, results in an appreciable value of r while the others predict a vanishingly small 
value. Future missions such as PLANCK may be capable of measuring this quantity with 
far greater precision. 

Summary 

We obtain a scalar spectral index Ug consistent with WMAP5 data by including fermion- 
dominated one-loop radiative corrections in a non-susy hybrid inflationary model. A com- 
pelling candidate for the origin of these corrections is a Yukawa coupling involving the infla- 
ton and right-handed neutrinos. A red-tilted spectral index is obtained for sub-Planckian 
values of the inflaton field, which was excluded to 2-cr by WMAP5 in the case of the tree 
level hybrid model. Furthermore, these corrections make accessible both non-hilltop and 
hilltop-type solutions in this model, and a red-tilted spectrum is achieved primarily for 
hilltop solutions. 
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